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We present a measurement of the time-dependent CP-violating asymmetries in B — » 
K*°j(K*° -» JfgTT ) decays based on 124 million T(4S) -> B-B decays collected with the BAR4R 
detector at the PEP-II asymmetric-energy B Factory at the Stanford Linear Accelerator Cen- 
ter. In a sample containing 105 ± 14 signal decays, we measure Sk*i = 0.25 ± 0.63 ± 0.14 and 
Ck*j = —0.57 ± 0.32 ± 0.09, where the first error is statistical and the second systematic. 

PACS numbers: 13.25.Hw, 13.25.-k, 14.40.Nd 



The recent data[lj from the B factory experiments 
have provided strong evidence that the quark mixing 
mechanism in the Standard Model (SM), encapsulated 
in the Cabibbo-Kobayashi-Maskawa (CKM) matrix^, is 
the dominant source of CP violation in the quark sec- 
tor. Nonetheless, decays which originate from radiative 
loop processes, such as b — > sj, may exhibit significant 
deviations from the SM due to new physics contribu- 
tions. In this letter we report the first measurement 
of time-dependent CP-violating (CPV) asymmetries in a 
b — > S7 process through the exclusive decay B° — > if*°7, 
where K* — > if°7r°[g. D. Atwood, M. Gronau and A. 
Soni were the first to point out that such a measure- 
ment probes the polarization of the photon jij, which 
is dominantly left-handed (right-handed) for b — ► 
{b — > s"f) in the SM, but is mixed in various new physics 
scenarios. The exclusive decays B° — > (K®ir )"fR and 
B° — > (K®ir°)jL are orthogonal transitions and are the 
dominant decays in the SM. Therefore the CPV asymme- 
try due to interference between decays with or without 
mixing is expected to be very small, ~ 2(to s /to;,) sin 2/3 
(/3 = &rg{-V cd V* b /V td V* b )). Any significant deviation 
would indicate phenomena beyond the SM. 

The B° — > if *°7 decays have been previously explored 
by the CLEO0, BAR4R@, and Belle collaborations 0, 
who reported measurements of branching fractions and 
the direct CP and isospin asymmetries. The measure- 
ments reported in this letter are based on 124 million 
T(45) -> BB decays collected in 1999-2003 at the PEP- 
II e + e~ collider at the Stanford Linear Accelerator Cen- 
ter with the BABAR detector, which is fully described in 
Ref. For the extraction of the time dependence of 
B° — > K*°j(K*° — > K®ir°) decays, we adopt an analysis 
approach that closely follows our recently published mea- 
surement of CPV asymmetries in the decay B° — > K®tt 
9J. There we established a technique of vertex recon- 
struction for B decay modes to final states containing a 
K® — > 7r + 7r~ decay and other neutral particles, but no 
primary charged particles at the B decay vertex. 

We search for B° -> K*°~/(K*° -> K°tt°) decays in 
hadronic events, which arc selected based on charged 
particle multiplicity and event topology. We reconstruct 
K° — ► 7r + 7r~ candidates from pairs of oppositely charged 
tracks, detected in the silicon vertex detector (SVT) 
and/or the central drift chamber (DCH). We require 
that these tracks originate from a vertex which is more 
than 0.3 cm from the primary vertex and that the re- 



sulting candidates have a 7t + 7t invariant mass between 
487 and 508MeV/c 2 . We form vr° -> 77 candidates 
from pairs of photon candidates in BABAR's electromag- 
netic calorimeter (EMC) which are not associated with 
any charged tracks, carry a minimum energy of 30McV, 
and possess the expected lateral shower shape. We re- 
quire that the 77 combination has an energy greater 
than 200 MeV and an invariant mass between 115 and 
155MeV/c 2 . We reconstruct candidate K* -> K°ir a de- 
cays from K®ir° combinations with invariant mass in the 
range 0.8 < M(if°7r°) < 1.0GeV/c 2 . For photons origi- 
nating from the B decay, we select clusters in the EMC 
which arc isolated by 25 cm from all other energy de- 
posits and are inconsistent with tt° — > 77 or 77 — > 77 
decays. 

We identify B° — > K*°j decays in if* 7 com- 
binations using two nearly independent kinematic 
variables: the energy-substituted mass to-es = 
\J (s/2 + p,j • ps) 2 /K 2 — v\ ano - tne ener gy difference 
AE = E* B - ^/s/2. Here (Ei,pi) and (E b ,Pb) arc the 
four-vectors of the initial e + e" system and the B candi- 
date, respectively, y/s is the center-of-mass energy, and 
the asterisk denotes th center-of-mass (CMS) frame. For 
signal decays, the tties distribution peaks near the B 
mass with a resolution of « 3.5MeV/c 2 and AE peaks 
near OMeV with a resolution of « 50 MeV. Both toes 
and AE exhibit a low-side tail from energy leakage in 
the EMC. For the study of CPV asymmetries, we con- 
sider candidates within 5.2 < toes < 5.3GeV/c 2 and 
\AE\ < 300 MeV, which includes the signal as well as 
a large "sideband" region for background estimation. 
When more than one candidate is found in an event, we 
select the combination with the ir° mass closest to the 
nominal 7r° value, and if ambiguity persists we select the 
combination with the K® mass closest to the nominal K ° 
value. 

The sample of candidate events selected by the above 
requirements contains significant background contribu- 
tions from continuum e + e~ — > qq (q — {u,d,s,c}), as 
well as random combinations from generic BB decays. 
We suppress both of these backgrounds by taking advan- 
tage of the expected angular distribution of the decay 
products of these processes. Angular momentum conser- 
vation restricts the K* meson in the B° — » K*°j decay 
to transversely polarized states, which leads to an angu- 
lar distribution of sin 2 8h for the decay products, where 
Oh is the angle between the K 2, and the B meson direc- 
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FIG. 1: Distribution of (a) jties and (b) Mk* for events en- 
hanced in signal decays. The dashed and solid curves rep- 
resent the background and signal-plus-background contribu- 
tions, respectively, as obtained from the maximum likelihood 
fit to the full data sample. The selection technique is de- 
scribed in the text. 



tions in the K* rest frame. Monte Carlo studies show 
that the background candidates peak near cos Or = — 1- 
We require cos Oh > —0.6, resulting in rejection of 68% of 
BB and 48% of continuum background candidates, while 
retaining 91% of the signal. 

We exploit topological variables to further suppress the 
continuum backgrounds, which in the CMS frame tend to 
retain the jet-like features of the qq fragmentation pro- 
cess, as opposed to spherical BB decays. In the CMS 
system we calculate the angle 9g between the spheric- 
ity axis of the B candidate and that of the remaining 
particles in the rest of the event (ROE). While |cos#g| 
is highly peaked near 1 for continuum background, it 
is nearly uniformly distributed for BB events. We re- 
quire | cos 0g | < 0.9, eliminating 58% of the continuum 
events. We also employ an event-shape Fisher discrim- 
inant in the maximum-likelihood fit (described below) 
from which we extract the CPV measurements. This 
variable is defined as T = 0.53 — O.6OL0 + 1.27L2, where 
L, 
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|p*|| cos 6* 1 3 , pi is the momentum of par- 



ticle i in the CMS system and 9* is the angle between p* 
and the sphericity axis of the B candidate. 

The above selections yield 1916 B° -> K*°j(K*° -> 
K^tt ) candidates. We extract our measurements from 
this sample using an unbinned maximum-likelihood fit 
to kinematic (toeSj AE, and K* mass), event shape (!F), 
flavor tag, and time structure variables (described be- 
low). As input to the fit, we parameterize the probability 
distribution functions (PDF) describing the observables 
of signal and BB background events using either more co- 
pious fully-reconstructed B decays in data or simulated 
samples. For the continuum background, we select the 
functional form of the PDFs describing each fit variable in 
data using the sideband regions of the other observables 
where the qq background dominates. We include these 
regions in the fitted sample and simultaneously extract 
the parameters of the background PDFs along with the 



CPV measurements. We fit 105 ± 14 signal and 19 ± 15 
other B decays in the selected sample. This signal yield is 
consistent with expectations from the previous measure- 
ments of the branching fractions 0, 0, Q- Figure Q dis- 
plays the toes and Mk* distributions for signal-enhanced 
sub-samples of these events, selected using the PDFs em- 
ployed in the fit (see below). 

For each B° — > K*°j candidate, we examine the re- 
maining tracks and neutral particles in the event to de- 
termine if the other B in the event, B tag , decayed as a B° 
or a B° (flavor tag). Time-dependent CPV asymmetries 
are determined by reconstructing the distribution of the 
proper decay time difference, At = tap — U ag - At the 
T(4S) resonance, the distribution of At follows 
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(S f sin (AtAm d ) - C f cos (AtAm d ))}, 



where the upper (lower) sign corresponds to B tag de- 
caying as B° (B°), t is the B° lifetime, Amd is the 
mixing frequency, and Sf and Cf are the magnitude of 
the mixing- induced and direct CPV asymmetries, respec- 
tively. As stated above, in the SM we expect 5if* 7 ~ 
2(m s /mb) sin 2/9 ss 0.05. We expect Ck*^ = —Ax»-y, the 
direct CP asymmetry measured in the self-tagging and 
more copious B a — > K*°^(K* — > K + ir~) decay. 

We use a neural network to determine the flavor, T, 
of the -Btag meson from kinematic and particle iden- 
tification informational • Each event is assigned to 
one of five mutually exclusive tagging categories, de- 
signed to combine flavor tags with similar performance 
and At resolution. We parameterize the performance of 
this algorithm in a data sample (i3fl av ) of fully recon- 
structed B° — > D(*'-K + /p + /a~i decays. The average 
effective tagging efficiency obtained from this sample is 



= S c e c s (l-2w c Y 



= 0.288±0.005, where e% and w c 



are 



the efficiency and mistag probabilities, respectively, for 
events tagged in category c. In each tagging category, we 
extract the fraction of events (e^g) and the asymmetry in 
the rate of B° and B° tags in the continuum background 
events in the fit to the data. 

We compute the proper time difference At from the 
known boost of the e + e~ system and the measured 
Az = zcp — z tagl the difference between the reconstructed 
decay vertex positions of the B° — > K* 0r y and B tag can- 
didate along the boost direction (z). A description of the 
inclusive reconstruction of the B tag vertex using tracks in 
ROE is given in [1(| . Replicating the vertexing technique 
developed for B° — > K^tt decays 0, we determine the 
decay point z C p for B° K*°j{K*° K°tt°) candi- 
dates from the intersection of the K® trajectory with the 
interaction region. This is accomplished by constrain- 
ing the B vertex to the interaction point (IP) in the 
plane transverse to the beam, which is determined in 
each run from the spatial distribution of vertices from 
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two-track events. We combine the uncertainty in the IP 
position, which follows from the size of the interaction 
region (about 200 /im horizontal and 4/im vertical), with 
the RMS of the transverse B flight length distribution 
(about 30 /im) to assign an uncertainty to the IP con- 
straint. 

Simulation studies indicate that B° — ► K*°~/(K* a — > 
K®ir°) decays exhibit properties which are characteristic 
of the IP vertexing technique, namely that the per-event 
estimate of the error on At, a At, reflects the expected 
dependence of the zqp resolution on the flight direc- 
tion and the number of SVT layers traversed by its decay 
daughters. Though the fit extracts CV» 7 from all flavor 
tagged signal decays, we only allow 68% of these events 
contribute to the measurement of Sx'-y- This subset con- 
sists of candidates which are composed of K® decays with 
at least one hit in the SVT on both tracks and pass the 
quality requirements of a At < 2.5 ps and |At| < 20 ps. 
For 66% of this subset, both tracks have hits in the inner 
three SVT layers, which results in a mean At resolution 
that is comparable to decays with the vertex directly re- 
constructed from charged particles originating at the B 
decay point ^(J- I n the remainder of the subset, the 
resolution is nearly two times worse. 

We obtain the PDF for the time-dependence of signal 
decays from the convolution of Eq. ^ with a resolution 
function lZ(5t = At — Ai tr ue, fAt)- The resolution func- 
tion is parameterized as the sum of a 'core' and a 'tail' 
Gaussian function, each with a width and mean propor- 
tional to the reconstructed a At, and a third Gaussian 
centered at zero with a fixed width of 8ps |ic|. Us- 
ing simulated data, we have verified that the parame- 
ters of K{6t,aAt) for B° -> K*°~f decays and the BB 
backgrounds are similar to those obtained from the Pflav 
sample, even though the distributions of a At differ con- 
siderably. Therefore, we extract these parameters from a 
fit to the Bflav sample. We find that the At distribution 
of continuum background candidates is well described by 
a delta function convoluted with a resolution function 
with the same functional form as used for signal events. 
We determine the parameters of the background function 
in the fit to the B° -► K*°~f(K* -> K°ir°) dataset. 

To extract the CPV asymmetries we maximize the log- 
arithm of the likelihood function 

L(bf,Cf,I\ h ,.t h ,e q? ,a) - (N s +N BB +N q9 )l x 

N qq f qq e q {P qq {x t ,yi;a)] x 
Ylie^/oAtlNsil- fs)e c s V' s (y i ;C f ) 

AW 1 - JbbYbb'P'bb^ + 

where the second (third) factor on the right-hand 
side is the contribution from events with (without) 



At information. The vectors afj and yi represent 
the time-structure and remaining observables, re- 
spectively, for event %. The PDFs Vh{ii,Vi) — 
P h {rn^i)Ph{AE i )P h {T i )P h (M K ^ i )P c h ^At l \aAt^T t ) 
wdV' h (yi) = P /l (m ESi )P /l (A J E; i )P /l (^)P /l (M if , )i )P ft Ci (r i ) 
are the products of the PDFs described above for hy- 
pothesis h of signal (S), BB background (BB), and 
continuum background (qq). Along with the CPV 
asymmetries Sf and Cf, the fit extracts the yields 
Ns, N B g, and N q q, the fractions of events with At 
information fs and f q q, and the parameters a which 
describe the background PDFs. We determine e B and 
f gB - in simulated generic BB decays. 

The fit to the data sample 
yields S K , 7 = 0.25 ± 0.63 ± 0.14 and 

C K , 7 = -0.57 ± 0.32 ± 0.09, where the uncertain- 
ties are statistical and systematic, respectively. The fit 
reports a correlation of 1% between these parameters. 
The result for Cr-* 7 is consistent with a fit that does not 
employ At information. Since the present measurements 
of Ak* 7 |S> 0] are consistent with zero, we also fit 
the data sample with Ck*j fixed to zero and obtain 
SW* 7 = 0.25 ±0.65 ±0.14. 

The event selection criteria employed to isolate signal- 
enhanced samples displayed in Figure ^ are based on 
a cut on the likelihood ratio R = Vs/(Vs + Vbb + 
Vqq) calculated without the displayed observable. The 
dashed and solid curves indicate background and signal- 
plus-background contributions, respectively, as obtained 
from the fit, but corrected for the selection efficiency 
of R. Figure [5] shows distributions of At for B°- 
and P>°-tagged events, and the asymmetry AK*~){At) = 
[Ngo — Ngo] I [Nb« + N B o] as a function of At, also for 
a signal-enhanced sample. 

We consider several sources of systematic uncertain- 
ties related to the level and possible asymmetry of the 
background contribution from generic BB decays. We 
estimate the impact of potential biases in the determina- 
tion of the BB background rate to lead to a systematic 
uncertainty of 0.04 (0.05) on Sk'j {Ck*j)- We estimate 
an uncertainty of 0.12 (0.03) due to potential CPV asym- 
metries in the BB backgrounds and 0.02 (0.06) due to 
possible asymmetries in the rate of B° versus B° tags 
in continuum backgrounds. We quantify possible sys- 
tematic effects due to the vertexing method in the same 
manner as Ref. 0] , estimating systematic uncertainties of 
0.04 (0.02) due to the choice of resolution function, 0.04 
(< 0.01) due to the vertexing technique, and 0.03 (0.01) 
due to possible misalignments of the SVT. Finally, we 
include a systematic uncertainty of 0.02 (0.02) due to 
tagging asymmetries in the signal and 0.02 (0.02) due to 
imperfect knowledge of the PDFs used in the fit. 

In summary, we have performed a measurement of the 
time-dependent CPV asymmetry Sk*i and the direct- 
CP violating asymmetry Ck*^ from B a — > K*°-f(K*° — > 
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FIG. 2: Distributions of At for events enhanced in signal 
decays with B tag tagged as (a) B° or (b) B°, and (c) the 
resulting asymmetry Ai(*-y(At). The dashed and solid curves 
represent the fitted background and signal-plus-background 
contributions, respectively, as obtained from the maximum 
likelihood fit. The raw asymmetry projection corresponds to 
approximately 38 signal and 19 background events. 



K®ir°) decays. Our measurement is consistent with the 
SM expectation of very small CPV asymmetries. 
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